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Remarks on the formation and decay of 
multidimensional shock waves 

V. G. Danilov* 


Abstract 

In this paper, we present a formula describing the formation and 
decay of shock wave type solutions in some special cases. 

In [I] and [2j, in the quadratic and general cases of convex nonlinearity, 
we consider the process of formation of shock waves for scalar conservation 
laws in the one-dimensional case. Recall that, in the construction suggested 
in these papers, the key role is played by the function U\{x) determined by 
the implicit equation 

f'(u 1 (x)) = -Kx + b, (1) 

where K > 0 and b are constants, and f(u) is the nonlinear (convex) density 
of the conservation law. 

In the present paper, we generalize this construction to the multidimen¬ 
sional case. 

The main point is to generalize Eq. (1). Recall that the function U\(x) 
in (1) describes both the shock wave formation and the decay of a nonstable 
step function (a rarefaction wave type solution). 

The problems of formation and decay of step functions are closely related 
to each other: the change t —> — t allows one to use solutions describing the 
step function formation to construct solutions describing the step function 
decay, and conversely. 
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This procedure is described in detail (in the framework of the technique 
used there) in [I; in the scalar quadratic case. 

We consider the equation 


du 

dt 


d 




Z— 1 


0, 


( 2 ) 


where f t (u) are smooth functions. 

To Eq. (2) there corresponds the system of differential equations (equa¬ 
tions of characteristics) 

Xi = fl(u), x\ t=0 =x 0 , (3) 

U = 0, u | t=0 = u 0 (x o). 

We have the following obvious assertion. 


Lemma 1 


J = det 


dx 

dx 0 


, d 2 fi(u 0 ) dup 

du 2 dxpi' 
i =i n 


( 4 ) 


The proof readily follows from the relation 


d 2 J 
dt 2 


0. 


A generalization of (1) is based on the fact that the derivative of the left- 
hand side of (1) with respect to x is exactly the expression under the sign of 
sum in (4) in the one-dimensional case. 

Namely, suppose that two smooth surfaces Ti and T 2 , Ti fl T 2 = 0, are 
given in a connected simply connected domain 

We assume that Tj are determined by the equations x* = y* (.s), seDe 
R"- 1 , i — 1,2, j — 1,... ,n. 

We also assume that in there exists a solution Ui(x) of the problem 


u l 


Ti 


E 


d 2 fi(u 0 ) 
du 2 


M 


dui 

dxi 


+ K 


0, 


= U = const , U\ | = const , 


( 5 ) 


where K = K(s) > 0 is an unknown function which we seek together 
with Ui(x). It is clear that problem (5) is the required generalization of (1). 
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Remark 1 The condition K = K(s) means that the function K{s) is con¬ 
stant on the characteristics corresponding to Eq. (5). 


Remark 2 The solvbility of (5) means that the vector held f"(uo) is not 
singular, i.e., 

in o, 

which is an analog of the convexity condition. 


In the multidimensional case, the function Ui(xi ,..., x n ), i.e., the solution 
of (5), will play the same role as the solution of Eq. (1) in the one-dimensional 
case. 

Problem (5) is equivalent to the following one: 


VU = X.'W. (6) 

«i | T , 0 = u, i = 1,.... n. 

We have 

Ui = U — K(s)t. 

Let t 0 (s ) be such that A"(r 0 (s) ,s) e r 2 , then 


dXj 

dr 

du\ 

dr 


/"M, 

-K, 


K(s) 


u-< 

h) (s) 


For given U and Mq, the condition that K(s ) is positive implies restrictions 
on the direction of motion along the trajectories determined by (6), and the 
fact that problem (5) has a solution means that Ti and T 2 are sections of the 
bundle determined by the trajectories of (6). 

To be definite, we assume that U > u\]. Then for K > 0 the motion along 
the trajectories of (6) must occur from Ti to T 2 with increasing r. Otherwise, 
problem (5) does not have solutions. 

Next, by O - we denote the domain lying “before” Ti, i.e., the domain 
entered by the trajectories of system (6) for r < 0. By 0 + we denote the 
domain lying “after” T 2 , i.e., the domain entered by the trajectories of sys¬ 
tem (6) for r > r 0 (s). 

By H ± we denote the characteristic functions of the domains and by 
H° we denote the characteristic function of the domain hi \ \ 
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( 7 ) 


We consider the equation 


du x - dfj{u) 

dt ^ dxi 


and set 

u\ t=0 = UH~+u° 0 H + + u 1 {x)H°. (8) 

Next, we must define the concrete geometry of the problem. For example, 
we can assume that T i and T 2 are closed and T 2 is located in the interior of 
Ti, or conversely. However, we shall not do this, but simply assume that we 
are interested in the solution of problem (7), (8) in the domain where it can 
be obtained from the initial condition by using the characteristics. 

Clearly, it follows from Lemma 1 and the choice of the function U\(x) in 
(8) that the wave turns over on the trajectories of the characteristic system 
corresponding to (7). 

More precisely, the trajectories of the characteristic system corresponding 
to (7), 

^ = /i(«i), Xi \ t=0 = x 0 , i = l,...,n, (9) 

such that the point (xio, • ■ •, x n o) = X(s,r) belongs to the trajectory of 
system (6) for some fixed s and 0 < r < 7o(s) intersect at t — t 0 (s ) = 1/K(s) 
at the point x = x*(s), x*(s) = x(t 0 (s), X(s,t)). 

Thus, for t > max s f 0 ( s )> the evolution of the initial condition (8) gives a 
shock type solution of the form 

u = U + H(S{x,t))(u° 0 -U), (10) 

where S G C °°, H(s) is the Heaviside function, and the set S(x,t) = 0 is the 
shock wave front. 

Moreover, at the point x at which the jump occurs, we have the inequality 

u + — < 0, (11) 

where u + is the limit value of the solution calculated along the trajectory of 
system (6) as x —> x; here x corresponds to the value r > r, f corresponds 
to x, and u _ is determined similarly. 

Inequality (11) can be treated as the stability condition for the jump of 
the solution to Eq. (7) in the multidimensional case. 
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Of course, we here must take into account the above assumption on the 
direction of motion along the trajectories of system (7). 

It is easy to see that the limit u + can also be calculated along the vector 
fu( u o), and the limit w_ along the vector /"([/). 

We agree to denote the limit of g(x,t ) as x —> x for fixed t along the 
vector X by 

(X) lim g(x,t), 

1C— 

and the limit of g(x,t ) as x —> x along the vector X but in the opposite 
direction by 

(X) lim g(x,t). 

X^*X 

Then the stability condition for the jump (11) can be written as 

(/IK)) dm u(x, t ) - (/1(C/)) lim u(x, t ) < 0. (12) 

x—>x X^X 

Definition 1 A piecewise constant solution of Eq. (7) of the form (10) is 
said to be absolutely nonstable if the following inequality holds at all points 
zx G {s(x, t ) = 0} for a fixed t: 

(/IK)) dm u{x, t) - (f" u (U)) lim u(x, t ) > 0. (13) 

rr— x—>x 

It follows from the above that an absolutely nonstable jump must turn 
into a solution of the form (8). This construction is completely similar to 
the one-dimensional case. More precisely, in the one-dimensional case, this 
is described for the case of quadratic nonlinearity in [TJ. The same also holds 
for the case of general convex nonlinearity in the one-dimensional case [2j. In 
the multidimensional case, in fact, the above assumptions reduce the problem 
to the one-dimensional problem along the trajectories of system (6). 

We note that system (6) can be easily integrated: 


ui = U - K{s)t, (14) 

X*{t,s) = xl(s) + - ft(U - Kt)), i = 

Now let Xq = y 1 (so) be an arbitrary point on Ti. By Xq = X(t 0 (so), sq) we 
denote the point of intersection of the trajectory of system (6) with T 2 - Let 
x(Xq ,t), i = 1,2, be solutions of system (9) such that 
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(15) 


Then, by (14), we have 

*(A' 0 \ t) - x(X 2 0 , t) = Xl - Xl + t(f\U) - /'(«»)) 

In the construction of multidimensional nonlinear waves, an important 
role is played by the level surface T* of the solution [X|. These surfaces 
are determine by the relations T* = {t = 0( x )}, where 0( x ) is the desired 
unknown function whose zero-level surface is assumed to be given (in our 
problem, these are the surfaces Ti and T 2 ). 

Clearly, if, for example, T° = Ti, then T* is the set of the endpoints of the 
trajectories of system (9) starting on W at time t. In this case, the function 
0i (x) is the time required for the trajectory starting at a point X 10 G Tx to 
come to the point X. Similarly, we determine 0 2 (:r) and T^. 

Let us consider the expression 


[0i (x) - 02 (a;)] 


rj 


= t - 0 2 (a;) 


rj 


where the restriction means that x is a point on T]. We have 
X = A" 10 + 0x/'(£/) = X 20 + 02 f( u o), 


where X 10 G Ti and X 20 G T 2 are some initial point of the trajectories (9). 
We have 

A0o - X 20 = 0 2 /'K) - ipid’(U). 

Next, we have A0 O = X 1 ( s ) for some s. We denote X 20 = A0(r o (s),s). 
By (15), we have 


(02 - 01 


f\u) = X 02 - x 20 + - f(u° 0 )](Kt - 1). 

L A 


Now let t = 1/A" and x G T[ fl T0 Then X 20 = A 02 . In general, we have 

= {Kt-l)K ~\f{U)J'{U) - /'(tig)) 
v t 


(02 -0l) 


wmw 

where (,) is the inner product in M n . Similarly, we define the quantity 


(02 


01 ) 


Tit 
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Everything said above is an analog of the Introduction in [2J. In the 
present text, we restrict ourselves to this and only note that we have prepared 
everything necessary to construct the multidimensional analog of the weak 
asymptotic solution given in (2j, which describes the formation of a shock 
wave. We shall present this in detail in the next paper. 

In conclusion, we formulate the solution of the problem concerning the 
decay of a nonstable step-function. 

Suppose that the initial data for Eq. 2 have the form (10) for t — 0, 
and the vector fields /'(mq) and f"(U), f"{u []) are transversal to the 

surface r 0 = {S(x, 0) = 0}. Next, we assume that inequality (13) holds at 
the points of r 0 . Then there exists a t > 0 such that for t G [0, t\, the solution 
of Eq. (2) with the initial condition 

has the form (8). 
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